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Morphological image analysis is applied to the time evolution of the probability distribution 
of a quantum particle moving in two and three-dimensional billiards. It is shown that the time- 
averaged Euler characteristic of the probability density provides a well defined quantity to distinguish 
- - - between classically integrable and non-integrable billiards. In three dimensions the time-averaged 

, mean breadth of the probability density may also be used for this purpose. 
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I. INTRODUCTION 



In the theory of quantum chaos, the connection between the chaotic behavior of classical systems and specific prop- 
erties of their quantum mechanical counterparts has received much attention. The most studied issues in this context 
CO ■ are the statistical properties of nuclear and atomic spectra [|l| and of two- (recently three) dimensional billiards [^|-D. 
In the mid-eighties, a conjecture (which became known as the BGS conjecture ||]) was made that the quantal energy 
level spacing of (ergodic) systems whose classical equivalent exhibit chaotic behavior obeys the Wigner Gaussian- 
. Orthogonal-Ensemble statistics, known from random matrix theory, whereas for classically integrable systems the 
. Poisson distribution applies. Exceptions to this conjecture are known |^,^, but it is generally believed that these form 
^ ■ a set of measure zero. 

Another way to make a connection between the chaotic behavior of classical systems and specific properties of their 
quantum mechanical counterparts is to look at the eigenstates of billiards. Eigenstates of a stadium may exhibit 
structure that corresponds to the classical orbits |9|- |ll|] . These structures, called scars, have been directly related 



. to unstable periodic orbits |]12[-p,4[. Visual inspection of movies that show the time evolution of the probability 
^ ■ distribution of a quantum particle moving in two and three-dimensional (2D and 3D) billiards of different shapes 
ir^ . strongly suggest that these sequences of images may contain enough information to distinguish between classically 

■ integrable and chaotic systems |l5| . 

\ To our knowledge the morphology of the patterns of time (in) dependent probability distributions has not been 
studied. In this paper we propose a method, based on concepts from integral geometry, to perform a morphological 
image analysis (MIA) of the probability distribution. We use MIA to characterize the solutions of the time-dependent 

■ Schrodinger equation for various 2D and 3D billiards. We demonstrate that indeed, one can define geometrical 
I and topological descriptors that are useful to distinguish between classically integrable and chaotic billiards. Our 

^ motivation to analyze time-dependent instead of stationary probability distributions is that the former usually provide 
information about many eigenstates simultaneously. By following the system over a sufficiently long period of time, one 
collects information on all those eigenstates that have non-negligible overlap with the initial state. This "parallelism" , 
absent in the stationary-state approach reduces the computational work. 

MIA amounts to the characterization of the geometry and topology of any pattern or series of patterns by means of 
the so-called Minkowski functionals, known from integral geometry |l7|-[l9|]. These functionals are related to familiar 
' quantities: In two (three) dimensions they correspond to the covered area, perimeter and connectivity (volume, 
surface area, integral mean curvature and connectivity) of the pattern. In integral geometry the calculation of the 
Minkowski functionals is relatively straightforward and requires little computational effort. MIA has proven to be 
a powerful tool in different fields where the description of patterns is important: Statistical physics, to describe the 
morphology of porous media and complex fluids; cosmology, to analyze the large scale distribution of matter in the 
universe; chemistry, to describe the morphology of patterns in reaction diffusion systems [|o|; seismology, to describe 
the spatial complexity of regional seismicity realizations [|l]; polymer science, to identify and quantify the morphology 
of mesoscale structures in block copolymers [^,^. Furthermore, in the semi-classical limit the sum of exponentials 
of eigenfrequencies of a billiard-shaped drum can be written in terms of the Minkowski functionals of the billiard [|4| , 
and the Weyl hypothesis states that the number of energy levels up to an energy E can be written in a similar manner 
(see e.g. 

In section 2 we briefly review the theory of Minkowski functionals and explain the technique of MIA in more detail. 
We then apply MIA to the motion of a quantum particle in various two and three-dimensional billiards and show that 
an appropriate quantity can be defined to distinguish classically integrable and chaotic billiards. Finally, we study 
the transition to chaos for a system whose degree of chaos is parametrized. 
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II. MORPHOLOGICAL IMAGE ANALYSIS 



At the heart of MIA hes the computation of Minkowski functionals of an arbitrary (thresholded) image. Therefore, 
in this section we wih give some mathematical background on Minkowski functionals (closely following |2^]) and sketch 
the outline of an algorithm to compute these functionals. 

We first define the Minkowski functionals for compact convex sets. The parallel set Kr of a compact convex set 
if at a distance r is the union of all closed balls of radius r, the centers of which are points in K. The well-known 
Euclidean measures of Kr^ such as the area and volume, can be written as polynomials in r. For example, the parallel 
volume of a square S and a cube C with edge lengths a are given by 

Vr{S) = c? ^^ar^nr^, (1) 
Vr{C) = + 6aV + 37rar^ + 47rrV3. (2) 

In general, the volume v^^^ in d dimensions of the parallel set Kj. at distance r of is given by the Steiner formula 

4^\K,) = Y^l'l\w'f\KY, (3) 
(/=0 ^ ' 

where wI'^\K) are the Minkowski functionals. For the first three dimensions it can be shown that [|l^ 

d = 1 : W^^\K) = 1{K), Wi'^ = 2xiK), (4) 
d = 2: wI,'\k) = A{K), l^f ) = ^U{K), = 7rx{K), (5) 

d = 3: wI,'\k) = V{K), l^f = ^S{K), W^'\k) = ^TrBiK), Ws{3){K) = yxW, (6) 

where 1{K), ^i^)^ U (K), ^i^)^ S{K)^ B{K) and x(^) denote respectively the length, area, perimeter, volume, surface 
area, mean breath and Euler characteristic or connectivity number of the convex set K. The mean breath B is closely 
related to the mean curvature H = {Ri + R2)/2RiR2 where Ri and R2 are the principal radii of curvature [ p6[ . 
The Euler characteristic x is related to the Gaussian curvature G = I/R1R2 and is equal to unity for a convex body 
| p6[ . The mean curvature and the Gaussian curvature are two useful measures of the curvature of a surface If 
the local surface geometry has an elliptic, Euclidean or hyperbolic shape respectively then the Gaussian curvature is 
respectively positive, zero or negative [2^] . 

Up to now, we have limited the discussion to convex sets. In order to deal with the more general sets one encounters 
in characterizing images, we consider the convex ring 7^, which is the class of all subsets A which can be expressed as 
finite unions of compact convex sets Ki 

I 

A = yi^,. (7) 

The Euler characteristic x is defined as an additive functional on 7^, so that for A e 71 

x{A) = X ( U ^0 = E^(^*) - E ^(^* riKj) + ... + {-iy+\{K^ n...nKi), (8) 

\i=l / i i<j 



and where for compact convex sets K 



The Euler characteristic describes A in a pure topological way, i.e. without reference to any kind of metric. It 
can be shown that this definition of x permits all Minkowski functionals on IZ to be written in terms of x ^i^nd the 
kinematical density jl^]. The functional x ^s defined in integral geometry is the same as the Euler characteristic 
defined in algebraic topology [|l^: For d = 2, x equals the number of connected components minus the number of 
holes and in three dimensions x is given by the number of connected components minus the number of tunnels plus 
the number of cavities. For example, for a solid cube x = for a hollow cube X = 2 and for a cube pierced by a 
tunnel X = 0- For multiply connected structures the Euler characteristic is negative. For complex structures it is 
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FIG. 1. Snapshot of the probabihty density at t = 319 (in dimensionless units, see text) in the case of a bilhard in the form 
of a quarter circle (left) and its thresholded counterpart using a threshold = 10% (right). The background is colored white, 
the objects black and the forbidden region for the quantum particle gray. 



often difficult to identify the number of connected components, tunnels and cavities. However, MIA directly yields x 
as will be explained below. 

In order to calculate the morphological properties of a particular image we first convert the image to a black-and- 
white picture using a threshold. Pixels (voxels) with intensity below the threshold are considered to be part of the 
background and the others are building up the objects in the 2D (3D) picture. According to integral geometry, the 
morphological properties of the various objects building up the black-and-white picture can be completely described 



in terms of Minkowski functionals |17|. In order to calculate the Minkowski functionals in an efficient way we consider 
each pixel (voxel) as the union of the disjoint collection of its interior, faces (for the 3D case only), open edges and 
vertices. The values of [/, V, B and x fo^" these single open structures can easily be calculated [^. By making 
use of the property of additivity of the Minkowski functionals and the fact that there is no overlap between open 
bodies on a lattice, we compute the Minkowski functionals of the whole image by starting from a complete white 
image, by calculating the change in the number of interiors, faces (for the 3D case only), open edges and vertices 
when one pixel (voxel) is added to the 2D (3D) image, and this until all pixels (voxels) building up the picture are 
added. The Minkowski functionals can then be computed from 

d = 2 : A = ns^ U = —4ns + 2ne, x = ~ + (10) 
d = 3 : V = ric^ S = —Gric + 2n/, 2B = Sric — 2n/ + ng, X = + n/ — rig + n^, (11) 

where {ric) denote the number of squares (cubes) in 2D (3D), nj counts the number of faces (in 3D only), rie and 
fly denote the number of edges and vertices respectively. From (10) and (11) it follows directly that the Minkowski 
functionals contain information about local 4 (8)-point correlation functions in 2D (3D). A more detailed description 
of the algorithm and an example of computer code to compute the Minkowski functionals is given in Ref. p7|]. 



III. METHOD 



The application of MIA to the time evolution of the probability distribution of a quantum particle moving in a 
billiard is straightforward. We solve the time-dependent Schrodinger equation for a particle moving in a billiard by 
a stable and accurate numerical method [|8|. For practical purposes the results obtained are exact. A time series of 
snapshots displaying the probability density can easily be extracted from this data. A collection of digital videos can 
be found in Each pattern of the time series is converted into a black-and-white picture by applying a threshold. 
The result of applying a threshold of 10% to a representative image of the probability density is shown in Fig]^. For 
each black-and-white image in the time series we compute the Minkowski functionals using the algorithm described 
above and analyze the behavior of the Minkowski functionals as a function of the control parameters such as threshold, 
initial conditions of the wave packet, and shape of the billiard. We use Gaussian wave packets as initial states and 
study a variety of classically integrable and non-integrable billiards (see Figs.||, ||,0). In our numerical work we express 
lengths in units of a fixed wavelength A and rescale energy (and time, setting h = 1) such that a wave packet with 
average momentum 27r/A has an average kinetic energy of one [^. 

As will become clear from the examples given below, to distinguish between classically integrable and chaotic 
systems it is expedient to define the time-averaged Minkowski functionals 

X{T,e) = j^ j\{t,9)dt ; B{T,9) = j^j^B{t,9)dt, (12) 
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FIG. 2. The shapes of the various two-dimensional billiards studied. A: square, B: equilateral triangle, C: quarter circle, D: 
semi circle, a: perturbed square, b: perturbed triangle, c: perturbed quarter circle, d: perturbed semi circle, e: semi cardiod, f: 
triangle with sides which ratio is irrational. The billiards are colored white and the forbidden regions for the quantum particle 
are colored black. In most calculations the linear size of the square was taken to be 13 A. 

and the quantities 



In Fig. 5 we present results that illustrate the typical behavior of the Minkowski functionals as a function of time, 
keeping the threshold fixed. The left panel of Fig. 5 shows the perimeter U (t, 0) as a function of time t for the billiards 
(A, a) using a threshold = 5%. Billiard A (solid line) shows a behavior which is manifestly different from that of 
billiard a (dashed line): For times t > 100 the solid line displays large fluctuations, in contrast to the dashed line. 
Other pairs of billiards (B,b) etc. show very similar, characteristic fluctuations in the perimeter as a function of time 
(results not shown). Essentially the same behavior is found for the Euler characteristic xi^^t)- example is given 
in the right panel of Fig. 5 for the case of billiards (B,b) and = 15%. Hence, in conclusion, for integrable billiards 
U{t^6) and x{^i^) show large fluctuations while for non-integrable billiards these large fluctuations do not occur. 

Since the Euler characteristic x{^-> ^) is a global measure of the curvatures of the objects in the black-and-white 
image and provides direct quantitative information about the topology of the probability distribution (see [|l8[, pp. 
34, 112 and 113) we will study the behavior of x as a function of t and 6 in more detail. For each choice of the 
threshold 6>, the cumulative time average of the Euler characteristic X(T, 6>) approaches a constant value which we 
denote by X{6). A representative example of the behavior of X(T^6) is shown in FigJ^. In Fig.^we plot the results 
for X(0)^ for each of the billiards of FigJ^. As 6 approaches zero, the thresholded black-and-white picture becomes 
completely filled with black pixels, hence \\me^QX{0) = A^/A, for any image. Obviously, for 6> ^ 0, X{0) does not 
contain any useful information and for reasons of clarity we therefore omit data for ^ < 1 in Fig. [^. From Fig.^ it 
is clear that the X{0) curves form two clusters, depending on whether the system is classically integrable or not. In 
fact, /i(x), can be used to classify billiards according to their classical (non)integr ability, as shown in Table 1. The 
maxima for the classically integrable systems all lie in the range [0.97, 1.10], whereas the maxima for the classically 
non-integrable systems lie in the range [1.29,1.38]. Further analysis indicates that the maximum of X{0) changes 
linearly with the energy of the wavepacket (results not shown). This scaling behavior could be trivially incorporated 
into the definition of X{0). Calculations (results not shown) for systems up to 19 A x 19 A suggest that the system 



X{0)= lim X(T,6>) ; B (i9) = lim B(T, 6>), 

T — ^oo T — >oo 

/i(x) = maxX(l9) ; /i(5) = maxB(l9), 

6 6 

where A = }? j A (A = A^/V) in two (three) dimensions. 



(13) 
(14) 



IV. RESULTS 2D 
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FIG. 3. The shapes of the various three-dimensional billiards studied. K: cube, L: sphere, M: cylinder, N: one eight sphere, 
k: perturbed cube, 1: perturbed sphere, m: perturbed cylinder, n: perturbed one eight sphere. In most calculations the sides 
of the bounding cube (not shown) were taken to be 9.5A. 

□ n 

a b 

FIG. 4. Two instances of the lemon-shaped billiard, with boundaries in the x — y plane defined by y{x) = ±(1 — 

1,1]. A:^ = l, a: ^ = 1.5. The billiards are colored white and the forbidden regions for the quantum particle are colored 
black. In most calculations the linear size of the square was taken to be 13 A. 

size-dependence of X{0) is very weak. Therefore, we conclude that the scaling properties /i(x) can be employed to test 
for the (non)integr ability of billiards. In contrast to what Fig.|^(left) might suggest, the cumulative time average and 
its infinite time limit of the perimeter U (and also the area A) do not have the same properties as X(T, 0) and X{0) 
(results not shown). Hence in two dimensions only /i(x) can be used to test for the (non)integr ability of billiards. 

V. RESULTS 3D 

For three dimensional billiards not much research has been done on the issue of characterizing classical chaos on 
basis of some quantum mechanical property. Very recently results of numerical work on the stability of classical 
trajectories [|2^, and on 3D Sinai billiards ^ has been reported, but much of the field is unexplored. 

In order to look for a quantum signature of chaos in 3D billiards (see Figjs]) we apply an analogous approach to the 
one described above for the 2D billiards. In three dimensions both the mean breadth B and the Euler characteristic 
X may be used to distinguish integrable and non-integrable billiards. From Figs.||and ^ it is clear that the X{6) and 
B(^) curves form two groups depending on whether the billiards are integrable or not. As was the case for the 2D 
billiards, we can again use the maximum values of these curves to classify the billiards according to their classical (non) 
integrability, as shown in Tables 2 and 3. For classically integrable systems /i(5) lies in the range [1.60,192] whereas 
for non-integrable systems jii{B) lies in the range [2.47,3.10], a significant difference. For the Euler characteristic, 
these ranges are [0.37-0.59] for integrable systems and [0.80,1.34] for non-integrable ones. Here the separation is less 
in magnitude but still present. 

VI. TRANSITION TO CHAOS 

An interesting case to confirm that /i(x) can be used to detect the transition from integrable to chaotic behavior is 
the lemon-shaped billiard B. Two instances of the lemon-shaped billiard are shown in Fig.0. The boundaries in the 



5 




FIG. 5. Left: Perimeter U(t, 0) (in pixels), normalized to the area A of the billiard (in pixels), as a function of time t using a 
threshold — 5%. Solid line: Billiard A. Dashed line: Billiard a. Right: Euler characteristic x{t^^)j normalized to the area of 
the billiard (in pixels), as a function of time t using a threshold 6 — 15%. Solid line: Billiard B. Dashed line: Billiard b. Time 
is measured in dimensionless units (see text). 
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FIG. 6. X{T,6) as a function of time T, using a threshold 6 — 15%. Solid line: Billiard B. Dashed line: Billiard b. 

X — y plane are defined hy y{x) = ±(1 — l^l"^), x e [—1, 1]. The integrability of the lemon-shaped billiard depends on 
the shape parameter S: For ^ = 1 and ^ = oo the billiard is integrable and otherwise it is not. If S slightly deviates 
from one, the billiard is classically partially chaotic (periodic orbits and chaotic trajectories coexist). Fig.|l^ shows 
/i(x) for the lemon-shaped billiard, for the shape parameter S varying between 0.5 and 2.0. The sharp peak centered 
around 6 = 1 shows that /i(x) is very sensitive for the presence of classically unstable trajectories. If there are any 
unstable trajectories in the classical system, the value of /i(x) foi" the quantum mechanical equivalent can be classified 
as chaotic. Thus, the curve in Fig. ^ clearly shows the transition from chaotic to integrable behavior and vice versa. 

VII. CONCLUSION 

We have described a method, based on concepts of integral geometry, to analyze the time evolution of the probability 
distribution of a quantum particle moving in two and three-dimensional billiards. We have demonstrated that the 
time- averaged Euler characteristic can be used to classify billiards as integrable or not. For three-dimensional billiards 
also the time-averaged mean breadth (or integral mean curvature) of the probability density can be employed for this 
purpose. For lemon-shaped billiards we have shown that the transition from integrable to chaotic behavior can be 
determined from the dependence of the time-averaged Euler characteristic on the parameter that controls the shape of 
the billiard. We believe it may also be of interest to apply the method described in this paper to the time-independent 
solutions. 
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FIG. 7. X(0) as function of the threshold (in percent), for all billiards shown in Fig. 



TABLE I. /i(x) and the corresponding for all billiards shown in Fig. 



Billiard 


e 


Kx) 


A 


10 


1.00 


B 


13 


1.10 


C 


13 


1.09 


D 


13 


0.97 


a 


15 


1.38 


b 


17 


1.38 
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16 


1.33 


d 


17 


1.29 


e 


18 


1.36 


f 


17 


1.34 
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FIG. 8. X(0) as a function of the threshold (in percent), for all billiards shown in Fig. 
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FIG. 9. B(^) as a function of the threshold (in percent), for all billiards shown in Fig. 
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FIG. 10. fi{x) function of the shape parameter S for the lemon- shaped biUiard. For (5 = 1 the biUiard is integrable. 
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